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Abstract 



In this paper we reconcile two contradictory statements about deep inelas- 
tic scattering (DIS) in manifestly covariant theories: (i) the scattering must 
be gauge invariant, even in the deep inelastic limit, and (ii) the pole term 
(which is not gauge invariant in a covariant theory) dominates the scattering 
amplitude in the deep inelastic limit. An "intermediate" answer is found to 
be true. We show that, at all energies, the gauge dependent part of the pole 
term cancels the gauge dependent part of the rescattering term, so that both 
the pole and rescattering terms can be separately redefined in a gauge invari- 
ant fashion. The resulting, redefined pole term is then shown to dominate 
the scattering in the deep inelastic limit. Details are worked out for a simple 
example in 1+1 dimensions. 

25.30.Fj, 24.85.+P, 13.60.Hb, 03.65.Ge, ll.lO.Kk 
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I. INTRODUCTION 



Deep inelastic scattering (DIS) is a very important tool for investigating the structure of 
hadrons IQ . It involves the interaction between an energetic lepton beam (electrons, muons, 
neutrinos) and the respective hadron. 

DIS has been described (i) in the framework of the Quark Parton Model (QPM) using the 
light-cone formalism and (ii) with the aid of covariant field theories. The QPM is based on 
the assumption that all quarks (before and after scattering) are on-shell . In this manner 
gauge invariance is built-in. The leading contribution is the pole (or Born) term and the 
rescattering term is negligable in the deep inelastic limit. In descriptions based on covariant 
field theory, however, the intermediate quarks (or nucleons in studies of the deuteron [Q) 
are off-shell, and hence the pole term cannot, by itself, be gauge invariant. In this case it 
would appear that the pole term cannot be the leading term in the deep inelastic limit, and 
this has been a problem for covariant theories for many years. 

To avoid this problem, de Forest replaced the current by its gauge invariant part. He 
presented no proof, and the "de Forest prescription" has always seemed adhoc. The problem 
continues to trouble calculations which rely on the use of off-shell particles and dominance 
of the pole term. In a recent paper Kelly outlines three alternatives. Assuming that the 
photon has four momentum = (i/, 0, 0, q), these three alternatives can be summarized as 
follows: 

• de Forest prescription: leave the time component of the current (the charge) un- 
changed, and replace the third component by J'^ — > {v/q)J^ . 

• Weyl prescription: leave the third component of the current unchanged, and replace 
the time component by {q/v)J^. 

• Landau prescription: replace the four current by + {J ■ q/Q'^) q^-, 

where we use the SLAG convention that = — > for electron scattering. Note that 
each of these descriptions gives a different result for the charge and components. In the 
generalized Breit frame, where q^ = (0, 0, 0, Q), the de Forest and the Landau prescriptions 
are identical. One of the major conclusions of this paper is a detailed theoretical argument 
justifying the Landau prescription. 

In addition to this, we present a simple toy model for the structure functions of a compos- 
ite system. The structure functions we obtain from this model exhibit scaling, are dominated 
in the deep inelastic limit by the "gauge invariant part" of the pole term, and give a simple, 
qualitative description of the distribution function of valence quarks. We present the results 
from two versions. In one, the particles are taken to be scalars, and the scalar "nucleon" is a 
bound state of a scalar "quark" with unit charge and a scalar diquark with no charge. In a 
second, more realistic model, the composite spin 1/2 nucleon is a bound state of a spin 1/2 
quark and a scalar diquark. Stimulated by the success of QCD in 1+1 dimensions 0, we 
carry out these model calculations in 1+1 dimensions, obtaining finite results without the 
need for model dependent form factors. The only parameters in the model are the masses 
of the quark and diquark, mi and m2 respectively. The 1+1 dimensional calculations do 
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FIG. 1. Elementary qd interaction. 



not necessarily give the right physics, but help in developing the necessary tools for a more 
complete treatment. 

The remainder of this paper is divided into five sections. First we present our simple 
bound state model for the nuclon (which could also be applied to the description of mesons). 
Then we study gauge invariance and show how to uniquely define the gauge invariant part 
of the pole term. In Sec. IV we show for scalar quarks that the gauge invariant part of 
the pole term does indeed dominate DIS. The model deep inelastic structure functions are 
calculated and compared with experiment in Sec. V, and conclusions are given in Sec. VI. 



We begin the discussion by constructing a dynamical model for the nucleon. In the first 
subsection we will assume that the scalar "nucleon" is a bound state of the two "fundamen- 
tal" scalar particles: a charged "quark" q of mass mi and an uncharged "diquark" d of mass 
m2. In the second subsection we will generalize to a charged spin 1/2 quark and a charged 
spin 1/2 nucleon. 



The two " fundamental" constituents will interact with each other via the coupling shown 
in Fig. 1: 



where A is the coupling strength, and the form factor is a function of the square of the 
diquark four-momentum. For generality, we keep these form factors for now, but they will 
be set to unity later when we carry out calculations in 1-1-1 dimensions. Hence the scattering 
amplitude is: 



II. BOUND STATE EQUATION FOR THE "NUCLEON 



A. Scalar "nucleons" as bound states of scalar constituents 



V = xfipDfip'i) 



(2.1) 



M = \f{pl)f{p'^) + {-if\'f{pl)f{p'^) h{s) + ... 



A/(pi)/(pf) 

l-\h{s) 



(2.2) 



where we have introduced the function h{s) defined by 
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FIG. 2. The bubble diagram. 



(2.3) 



{2nY (ml — /c^ — ze)(mf — (p — /c)^ 



with p = P1+P2 and s = p^. This "bubble" integral is represented graphically in Fig. H. The 
mass of the bound state must be below threshold, and for simplicity we assume m2 — mi < 
M < mi + m2. 

In this simple separable model, the bound state is described by the "nucleon"g(i vertex 
function 



where A/" is a constant to be determined from the normalization of the wave function. This 
is shown in Fig. |^. Denoting the dressed propagator of the bound nucleon by Ab(p), the 
dressed qd scattering amplitude is given exactly in this simple model by 



V{p,) = Mfipl) , 



(2.4) 



M 



V{p2)^b{pW2) 



(2.5) 



where 



M^Ab{p) 



A 



(2.6) 



1- A6(s) ■ 



Eq. ( p3D is illustrated in Fig. |. 




FIG. 3. The "nucleon" vertex. 
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FIG. 4. The dressed bound state pole contribution to the qd scattering amplitude. 

Since the propagator should have a pole at s = M^, it is useful to expand the bubble 
h{s) in this vicinity, so the denominator of the RHS of Eq. (|2.6| ) is 



1 - Xh{M^) + A(M2 - s) h\M^) + ■■■ 

and this must vanish at s = M^. Consequently A6(M^) = 1 is a necessary condition of the 
bound state, and near the pole Eq. ( p.6|) becomes 



Knowing that, in the vicinity of s = M^, the nucleon propagator should be 
allows us to find the normalization of the Nqd vertex function: 



Note that the existance of the bound state implies that b'{M'^) < 0, which will be shown 
below. 

In 1+1 D the bubble integral will converge without form factors, and we can assume, for 
simplicity, that /(P2) = 1- The bubble integral ( p.3| ) is then model independent, and easily 
calculated. If it is parametrized using Feynman variables, one can either compute the loop 
below treshold directly in Minkowski space (doing the momentum integration first and then 
integrating over the Feynman variable) or compute it in Euclidian space and Wick rotate 
back to the Minkowski space. Introducing the triangle function 

A{p'^,ml,ml) = {p^ - ml- m\f - ^m\m\ , (2.10) 

which is negative if (m2 — rriiY < < ("^2 + ^1)^, the bubble becomes: 



o{s) = —- j^^^=^^^= arctan 

2^ J—A{s, mf, ml) 



-A(s, ml, ml) 

2 2 
— mf — m2 



(2.11) 



This is real, and if the mass of the bound state is to be M, then A must be 
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-27r-y-A(M2, mf,m 



arctan 



-A(M2,mf,mi) 
— mf — ml 



(2.12) 



This is real (so the Lagrangian is hermitian, as it should be) and negative (as expected from 
the fact that we have a bound state, implying the interaction is attractive). In fact, the 
whole 1+1 D theory can be derived from the Lagrangian 



L = -{d^d){d,d) - -mid' + {d^q){d,q)^ - mlqq^ 
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(2.13) 



where q is the "quark" and d the "diquark" fields. From this Lagrangian, or from Eq. fp.ll 
one sees that A must be negative in order for a bound state to exist. 
Above threshold, however, the loop h{s) also has an imaginary part 



96(s) 



-1 



2JA(s,mf,mi) 



(2.14) 



which can be calculated in two ways. The first method is a computation based on Cutkosky's 
rules, which are easily handled in this simple case. The second follows directly from the 
Feynman parametrization of Eq. (|2.3|) . Integrating over the internal momentum gives 



h{s) = - 



dx 



x] 



An Jo rriiX + 7712(1 — x) — s x{l 
Introducing the variable s' = m\/x + m\/{l — x) transforms ( |2.15| ) into: 



(2.15) 



b{s) 



1 



ds' 



TT J(mi+m2)2 S' — S — 16 



2JA(s',m?,mi) 



(2.16) 



This expression yields a real value for b{s) if s is below threshold because A{s' ,mf,m2) is 
positive above threshold. This relation is a dispersion relation, and comparing it with the 
general form 



m = i / 

TT Jti 



ds'^fis') 



troshold S' — S — le 



(2.17) 



one concludes that the imaginary part of the bubble is correctly given in Eq. ( |2.14| ). As a 
byproduct, we have also discovered the correct phase space factor for use with dispersion 
relations in 1+1 D. 

The real part of b{s) above treshold can also be found in two ways: (i) using Eq. ( p.lTj) , 
or (ii) by analytically continuing Eq. ( |2.11D , a procedure which also provides a third method 
for computing the bubble's imaginary part. As a result of these calculations, the real part 
above treshold is found to be: 



In 



47rw A(s, ml, ml 



— s{ml + ml) + sJA{s, mj, ml) 



— s(mf + ml) — sJA{s, mf, ml) 



(2.18) 
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At large s, Qb{s) ~ s~^, ^b{s) ~ s~^ln(s/m|); so b{s) — > as s — » oo. This important 
result will be used below. 

Before turning to the spin 1/2 case, we use the representation ( p.l5| ) to write the nor- 
malization condition (|2.9| ) in the following form: 

, x(l- x) , , 

1 = / dx ^ ^. 2.19 

47r 7o [mjx + ml{l - x) - M'^ x{l - x)Y 

This relation will be used again in the following sections. 



B. Generalization to spin 1/2 nucleons and quarks 

In this subsection we will repeat our former calculations for the more realistic case of 
charged spin 1/2 nucleons and quarks (retaining, however, the uncharged scalar diquark). 
The quark-diquark contact interaction is assumed to be a scalar in the Dirac space of the 
quarks, with the momentum space dependence previously given in Eq. ( p.l|) . 

In this case the bubble in d dimensions becomes: 

b{p) = I / . ... , . , , ' J — , 2.20 

^ ' J (27r)'^ (mi - P - ze) (mi -(/»-/;)- ze) ^ ' 

or, after the integration: 

b{p) = Aip") + B{p^) A (2.21) 

In 1+1 dimensions with the form factors set to unity, the factors A[p'^) and B{p'^) become: 
, , 9, mi dx 

Bip'^) = — ^ / — 7 ^'^ ^ , , (2.22) 

47r Jo m\x + m2(l — x) — p^x\ \ — x) 

where x is the Feynman parameter used previously in Eq. ( p^.l5| ). In the presence of spin, 
Eq. ( p.2| ) is modified: 

^ A (1 - \A{p^) + \B{p^) j)) 

(1 - AA(p2))2 _ X2B^{p^)p^) ■ ^ ^ ^ 

Expanding the denominator near the nucleon pole gives 

(1 - AA(/))2 _ X^B\p^)p^ = (1 - XAof - 2XA'{1 - AAo)(p' - M^) 

-X\B^o + 2BoB\p^ - M^))/ + . . . , (2.24) 

where Aq = A{p'^), Bq = BijP'), the primes denote derivatives with respect of p^, and all 
functions are evaluated at p^ = M^. The denominator vanishes at p"^ = if 

1 - AAo = MXBo. (2.25) 



4iT Jo mix + m^l — x) — p'^x{l — x) 



2s ^ (1 ^) dx 
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FIG. 5. The pole term. 

This is the bound state condition. After imposing the bound state condition, in the vicinity 
of the poles the scattering ampHtude becomes: 

M ^ (2 26) 

The normahzation factor M is then: 

J\f-^ = - [2M{A' + MB') + Bo] . (2.27) 



Using the relations ( |2.22|) , and the identity 



/■I , m?x^ — mo(l — x)^ 
0= / dx— \ 2. (2.28) 

Jo [rriiX + — x) — M'^x{l — x)] 



we may rewrite the renormalization condition ( |2.27|) in the following form: 



47r io [mjx + ml{l - x) - M^x{l - x)f 
This relation will be needed later to verify baryon number conservation. 



III. DEEP INELASTIC SCATTERING AND GAUGE INVARIANCE 

We now add electromagnetic interactions to the model, and study the implications of the 
requirement that the e.m. interaction be gauge invariant. We begin by restoring the form 
factors; later we will move to 1+1 dimensions and set the form factors to unity. As in the 
previous section, we will first carry out the discussion for a scalar quark and nucleon, and 
later extend it to a spin 1/2 quark and nucleon. 

A. The scalar model 

Start with the pole term, shown in Fig. ^, which gives 

eJJ^ = t{-tYenp,,p[)Aip[)Tip,) , (3.1) 
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FIG. 6. The rescattering diagram. 



where p'^ = Pi — q, the propagator of the scalar (charged) "quark" is 

A(pi)= 2 \ ■ ^ (3-2) 
mf — pf — le 

and j'^{pi,p'i) is the one-body current of the quark. (Note that our definition of the current 
JT" does not include the charge.) We assume that the diquark has no electromagnetic inter- 
action, so there are no diagrams with the photon coupling to the diquark. We require that 
the quark current satisfy the Ward-Takahashi (WT) identity: 

g%(Pi,K) = A(K)-A(pi). (3.3) 

The current of a bare, spinless particle 

j''iPup[) = iPi+p'ir (3.4) 

satisfies this identity. 

Using the WT identity, and noting that A^^(pi) = because pi = ml, we see that the 
pole term by itself is not gauge invariant 

q,Ja^ = -q,npi,p[)^{p[)np2) 

= -'A-\p[)-A-\p,)] A{p[)T{p,) 

= -T{p2)^0. (3.5) 

Hence the pole term by itself cannot be a good approximation for the scattering amplitude. 

To obtain a gauge invariant result, it is necessary to add the rescattering term shown in 
Fig. ^. This term is 

^6. A/(pi) 



d'k f{k^)np'-k,p-k)T{k) 

(27r)'^ {ml — {p' — ky — ie){ml — {p — ky — ie){m2 — k"^ — ie) ' 

where p' = p + q and s' = p'"^ is the square of the momentum of the final state after it has 
absorbed the virtual (space- like) photon. Note that p'^ = s = but that p'^ = s' > M^. 
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This diagram is also not gauge invariant. Using Eq. ( |2.4| ) and the WT identity gives 



1 - Xb{s') 



Xh{s') - Xh{M^) =Mf{pl) = T{p, 



(3.7) 



where the bound state condition A6(M^) = 1 was used in the second step. This is not gauge 
invariant, but if we add the rescattering term to the pole term the total amplitude is 



(3.8) 



Stated in a different way, the gauge dependent parts of the pole term and the rescattering 
term cancel. 

We now face the question posed in the introduction: How can the pole term dominate 
the scattering in the deep inelastic limit if it is not gauge invariant? To answer this question 
we break the pole term into a gauge invariant part and a gauge dependent part as follows: 



where 



(2pi - g)^r(p2) 
mf — {pi — qY 



' 2(g-p0 



2 



{pi - qf 



-q' 1 



q' 



{pi - qf - ml 



Note that g^/f = 0, and that the gauge dependent part of the pole term reduces to 



q 



^r(p2). 



q^ 



Now the rescattering term can also be similarly decomposed. Noting that 
f^p> -k,p-k) = If (p' -k) + I^ip' - k) 
the gauge dependent part of the rescattering term becomes 



(f_ \T{p2) 

'g2 [1-A6(s')] 

df'k P{k^){{p-k + qf-{p-kf] 



(27r)°' {ml — A;^)(mf — {p — A;)^)(mf — {p — k + qy 



q_ 

q' 



g2 



r(P2 



(3.9) 

(3.10) 
(3.11) 



(3.12) 



(3.13) 



(3.14) 



which exactly cancels the gauge dependent part of the pole term! 

We have shown that, if we drop the parts from each terms, the Born pole term and the 
rescattering term are redefined so that each is separately gauge invariant. As it turns out, 
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the redefined pole term is identical to the Landau prescription defined in the introduction. 
We have just justified this prescription. 

It is clear that the same decomposition works in 1 + 1 dimensions, where the form factor 
/ = 1. In this case the diagrams come from from the Lagrangian (|2.13|) , modified to include 
electromagnetic interactions. The new Lagrangian is 

^ = lid^d){d,d) - \mlS + {D,q){D>^q)^ - mlqq^ - ^d'qq^ , (3.15) 

where the fields are the same as in the preceding section, but the gradient of the charged 
field has been replaced by the covariant derivative 

D^" = d^" - leA^" . (3.16) 



B. Extension to sin 1/2 quarks and nucleons 



We how extend the above results to the spin 1/2 model previously introduced. The Born 
term for a structureless quark (with the form factor set equal to unity) is 



eJa' = -ef{p,,p[)Sip[)T{p2,p) 
where now j'^ = and the quark propagator is 

S{Pi) = — -• 

mi— p' — te 
The vertex function for the spin 1/2 case is 

r(p,A;v) =AruN{p,XN) 



(3.17) 



(3.18) 



(3.19) 



where Atv is the helicity of the nucleon and Un{p, A at) are the helicity spinor of the nucleon. 
Using the Dirac equation for the final quark, the pole term can be written 



J' 



Af 



2 

mf 



N 



(3.20) 



where Ai the helicity of the outgoing quark and u{pi,\i) the hehcity spinor of the quark. 
Now the current is spin dependent. 

As before, we decompose the pole term into a gauge invariant and gauge noninvariant 
part: 



m\ — {pi — qY 



u{pi, Ai 



2pi ■ qqf" 



Un{p, Atv) 



Jag = -— A/'m(pi, Xi)Un{p, Xn) = m(Pi, Ai)r(p, A 



q- 



N 



(3.21) 



Note that, except for the addition of the spinors, the gauge dependent part identical to the 
scalar case, and the gauge independent part is again the Landau prescription. 
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We now look at the spinor rescattering diagram. Introducing the decomposition 



" - ^ - ^ {S-\v -k + q)- S-\p - k)} , (3,22) 



q^ q 

the gauge dependent part of the rescattering term becomes 
Jl^ = ^ X u{pu Xi) [1 - Xb{p + q)]-' 



g2 



uipu Ai) [1 - Xbip + q)]-' {a b{p) -Xbip + q)} Tip, Xn) 



= ^u{pi,Xi)T{p,Xn), (3.23) 

where we used the bound state condition ( |2.25| ) and /^r(p, Aat) = MT{p, Xn)- Note that 
this term exactly cancels the gauge dependent part of the pole term! 

This argument shows that this technique is quite general, and should work in every case. 
The pole term is divided into a term which is gauge invariant and one proportional to q'^ by 
adding and subtracting a (unique) term proportional to q^. Then it may be shown that this 
non gauge invariant term is cancelled when the rescattering term is similarily separated into 
(unique) gauge invariant and non guage invariant parts. Knowing that this works means 
that the non gauge invariant part of the pole term can be discarded without checking that it 
is cancelled by an identical piece from the rescattering term. Our argument can be extended 
to currents with form factors if we use the technique developed by Gross and Riska to 
define the off-shell currents. 



IV. RESCATTERING AND THE TRIANGLE DIAGRAM 

In this section we consider the rescattering term for the scalar model in more detail. We 
compute it exactly for the case g = and confirm that the charge of the bound state is 
properly normalized. Then we study the behavior of the elastic form factor and the inelastic 
rescattering term at large g^. We will compare the gauge invariant part of the rescattering 
term with the gauge invariant part of the pole term and show that the pole term dominates 
at high g^. In order to achieve convergence in a model independent way, we return to 1+1 
dimension. 

The key to these calculations is the triangle diagram. The gauge invariant part of this 
diagram is 

fNr>' r>\ = -- [ ^ Q^P-ky-Q^'Q-ip-k) 

q^J (27r)2 (mi - P)(m2 - (p + g - A;)2)(mf - (p - k)^) ' ^ ' ^ 
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FIG. 7. The hadronic part of the elastic electron-nucleon scattering process, proportional to 
the deuteron form factor. 



= U - /" , (4.2) 



Feynman parametrizing the integral and doing the momentum integration gives 

' q- p\ 1 /"i-^ , I - X - y 

— - — — ax ay — 

J 8n Jo Jo D^{x,y,p,q) 

where the denominator is 

D{x, y, p, q) = {px + gx + py^ 

+x {m\ — {p + qY — ^2) + y ij^i ~ ~ "'^2) + "^2 • (4-3) 

It is convenient to express the integral in terms of p' = p + q and p, and to change to the 
variables 

z = X + y 

u = ^ , (4.4) 
which maps the region of integration into 

0< ;z < 1 

-l<u<l. (4.5) 
The Jacobian of this transformation is z/2, and the integral then becomes 

a'^ Sn Jo J -I d'^iz^u^v'.v) 



where 



d{z^ u,p',p) = m\{l — z) + 'm\z — \{p'^ + p^)z{l — z) 

-lip" - p')uz{l -z)- \q'z\l - u") . (4.7) 

We will use the general results ( [4.6| ) and ( |4.7] ) in the following discussion. 
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A. Charge Normalization 



First consider the hadronic part of the elastic electron-nucleon scattering process, shown 
in Fig. 1^. This diagram is gauge invariant by itself, and hence the deuteron form factor is 
given by 



ej^ = e F{q^) {p' + pY = eU^ tf'^p', p) 



(4.8) 



where, for elastic scattering, the triangle diagram is evaluated with the constraints p'^ 
{p + qf = M'^ = p^. Substituting Eq. (|]6]) into Eq. gives 



' Svr io 



du 



z(l-z) 



a 



(z) - Wz^{l-u^) 



2 ' 



where 

a{z) = ml{l - z) + m\z - M^z{\ - z) . 
If the charge is conserved, the form factor must satisfy the condition 

z(l-z) 



F(0) = 1 = — / dz du 

Stt jo J-1 



d 



9M2 47r 



— f dz 



[mi(l -z)+ m\z - M2z(l - z)\ 
1 



[mi(l - + m\z - M^z{l - z)] 



Comparing this with Eq. ( p.l5|) gives the relation 

1 = -Af' 



.^db{M^) 



9M2 



(4.9) 



(4.10) 



(4.11) 



(4.12) 



which is in agreement with the Eq. ( p.9|) and proves that charge is conserved. 



B. Behavior of the Form Factor at Large Momentum Transfers 

In preparation for the discussion of rescattering in the deep inelastic limit, and in order 
to study the properties of the 1+1 D model, we examine the behavior of the form factor in 
the high limit. This is obtained from 



d 



d 



— I dz 



du 



1 



a(z) + iQV(l-M2) 



1 dz 



8nQ Jo 



9M2 [87rg7o z /a{z) + \Q^z^ 



z)dz 



log 



a{z) + \QH^ + \Qz 
a{z) + \Q^z^ - \Qz, 



Ci{z) ^a{z) + \QH^ 



(4.13) 
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where we have introduced = —q^ > used in electron scattering. At large Q, the last 
integral is dominated by values of z near zero, and is well approximated by 



8nQ Jo c^(o) y^a(O) + \Q^z^' 

W r^^i^ 1 



log 



47rm2(y \m2/ 

We see that F{Q'^) falls off as Q approaches infinity, as it should for a composite particle. 



C. Implications for protons, neutrons, and mesons 

These results can be applied to realistic cases. Denoting the charge of the quark by Cq 
and the mass of the diquark by m2 = rrid the total charge is simply J2i Cj, as required by 
charge conservation. At high however, we obtain 

eW)^i:e.n,W')-^i:;|>o.(|-) . (4.15) 



Note that the sign of Fp{Q'^) is positive for the proton, but that it is identically zero (or, 
more correctly, falls off faster than for the neutron, unless the mass of the ud diquark 
differs from the mass of the dd diquark. In the latter case. 




yet been determined experimentally, and continues to be a subject of great interest. 

Note that Eq. ( [4.15|) also predicts that the electromagnetic form factors of spin qq 
mesons which have the same flavor for the quark and antiquark (such as ss, cc) , is identically 
zero. 



D. Rescattering in the deep inelastic limit 

In Sec. Ill the rescattering term was treated briefly. We divided the triangle into a gauge 
dependent and gauge independent part. The former was calculated exactly, while the latter 
is proportional to the triangle diagram given in Eq. (^.6| ): 
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FIG. 8. The two virtual particles in the final state can both be on shell (symbolized by the 
x's) if W > mi + 1712- 



hp'+pr-'-^^^^]nQ^w^), (4.17) 



[1 - Xb{W^)] ['^ ' q 
where W'^ = s' = p'"^ and the scalar function T{Q'^ , W"^) is 

— [ dz [ du ^ 

^TT Jo J -I [aw{z) + I3{z,u) - it] 

where 



T{Q\W^) = — I dz I du ^ \ (4.18) 



aw{,z) = m\{l — z) + m\z — W'^z{l — z) 

f3{z, u) = |(1V' - M^)z{l -z){l-u) + iQV(l - u^) . (4.19) 

For electron scattering, > and W'^ — > 0, so /3 > 0. However, as for the bubble 
diagram, aw has a zero ifW> mi + m2. This means that the denominator of Eq. ([4.18|) also 
has a zero, and that T has an imaginary part. From the principles of dispersion theory we 
know that this is due to the existance of physical scattering in the final state, as illustrated 
in Fig. |. 

To obtain a tractible form for the integral, replace the u integration by (3, which gives 



TiQ',W') = — dz d(3 I ^ -, (4.20) 

where 

7(z, = \[{W' - M')z{l -z) + Q'z']' - Q'z'p , (4.21) 
We now replace by 

W^ = {p + qf = + Q2 g _ 1 j (4.22) 
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where x is the Bjorken scahng variable 



^ = = . (4-23) 

2p-q 2Mz/' ^ ' 

and take the deep inelastic limit where Q"^ ^ oo with x held constant. The integral is 
dominated by values of z near zero, and hence may be approximated by 

T(Q2, VV) ~ — — ^ dz dp 



2ttQ^{1-x)Jo Jo [ml-W^z + l3-ie] 

X , f 1 l\ 

dz 



2nQ'^{l — x) Jo \[m2 — W^z — ie] ml^ 

X f inx 2 \ 

27iQ^{l-x) \Q^{l-x) " ' 



(4.24) 



where we have kept the leading contribution from the imaginary part even though it is 
negligible compared to the real part. 

In conclusion, one can state that the gauge invariant part of the rescattering term falls 
like l/Q^. Noting that the bubble also falls as l/Q^, the full rescattering term in the deep 
inelastic limit approaches 

This is to be compared with the deep inelastic limit of the Born term, which is 

= -2 (rf - ^±^A , f „ , (4.26) 

\ J mf- {pi - qy 

and does not fall with Q^. Hence, in the deep inelastic limit the gauge invariant part of the 
rescattering term is negligible in comparison to the gauge invariant part of the pole term. 

In the next section we will calculate the cross section for deep inelastic scattering, and 
find the structure functions predicted by this simple model. 



V. THE NUCLEON STRUCTURE FUNCTIONS 

In this section we calculate the cross section and structure functions for the deep inelastic 
scattering (DIS) of unpolarized electrons from the composite "nucleon" A*", 

e + N^e' + X, (5.1) 

where X is the undetected final hadronic state. This scattering process is illustrated in 
Fig. ^. The nucleon is composite, as described in Sec. II. For completeness, the well known 
kinematics and cross section formula [0,0 are first reviewed in the next two subsections. 
The remaining subsections give results for the 1+1 D toy models presented in this paper. 
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FIG. 9. The leading contribution to DIS from a hadron. 



A. The Kinematics of DIS 



The Bjorken variable, x, was defined in Eq. ( f4.23| ) and the invariant mass of the final 
hadronic state, W, in Eq. ( [4.22|) . In terms of these quantities, the energy transferred by the 
photon to the hadron in the lab system is 

gO = Zy = -9L (5.2) 

^ 2Mx ^ ' 

and the magnitude of the spatial component of g, which will be chosen to be in the z 
direction, is: 



|q| = 



In the CM of the outgoing hadronic pair, the invariant mass of the final state is 



(5.3) 



(5.4) 



where the prime denotes that the variable is evaluated in the CM system. The magnitude of 
the three-momentum of the "quark" (and "diquark") in the CM frame, and in the Bjorken 
limit, is 



,2 A{W^mlml) Q\l-x) 



Pi 



4iy2 



4x 



1 + 



x{M^ - 2ml - 2mi) 
Q'{l-x) 



(5.5) 



where the non-leading term will be needed below. The components of q in the CM frame 
can be found by a Lorentz boost from the lab: 



f<lo\ 




(\ 


+ ql 








-(Iz \ 




fqo\ 





1 







1 
















~ w 










1 










U'J 




[ 










JW^ + ql ) 







(5.6) 
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Hence the CM components are: 



2Mx 

^. = Q 



+ Q2 ( - - 1 



-1/2 



1/2 



+ — + 



2Mx 



2Mx 



1 



X 



-1/2 



X 



2 , , Q' 



1/2 



In the deep inelastic hmit these become 

Q{l-2x) 



2Mx 



(5.7) 



2^x{l-x) 

Q 

2Jx{l-x) V 



1 - 



2Q2(i 
M2x(1 - 2x)2 

2g2(i-x) 



(5.8) 



where the non- leading terms will be needed later. Using the same boost, the energy of the 
initial hadron in the CM of the outgoing hadronic system is 



M 



M2 



1 



X 



.-1/2 



' g2 Q 



1/2 



X 



4M2x2 



Q 



^ M'^x[l-A{l-x)] 
2^x{l-x) V 2Q2(i -x) 

as expected. We now turn to the calculation of the cross section. 





' M2' 


/ 

= <lz 





(5.9) 



B. Cross section and structure functions 

The scattering amplitude for the DIS process is 



e2 



Mf. = —^3e,{k',k)J^, (5.10) 

where is the hadronic current (kept general for now) and jefj,{k', k) is the current of the 
electron 

j,^ik',k)=u{k's'hXk,s). (5.11) 

The spins of the electron are denoted by s and s', and the momenta are labeled as in Fig. |^. 
In this notation the unpolarized differential cross section in the lab is 



/ dJdPidP2 (27r)^5(fc' + pi+p2-fc-P)|-M/^r 

Ja {27iy32MEE'EiE2^ ^ ^ ^/^i^i^i j\ m 

= A (27r)38EE'M ^'^'^ ' ^^"^^^ 
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where |A^/jp contains a sum over the spins of the final particles and an average over the 
spins of the initial particles, Ei and E2 are the energies of the outgoing quark and diquark, E 
and E' are the initial and final electron energies, and i'^" and W^" are the lepton and hadron 
current tensors, defined below. The hadronic tensor contains all of the the physical infor- 
mation about the hadron-photon interaction and includes the integration over the outgoing 
hadrons and the normalization factors associated with the hadronic wave functions. 

If the electron mass can be neglected the outgoing electron differential can be written 
d^k' = E''^dE'dQ' , which gives the following result for the inelastic cross section: 



g4 [e i AnM 



dQ' dE' 

If the electron mass is neglected the lepton current tensor is 



(5.13) 



J2 ^luJet, = 7,Tr (7^ /;'7^ /;) = 2 (k'^K + k^k'^ + - q^g^^ 



spins 



Gauge invariance of the hadronic current implies 



0. 



(5.14) 



(5.15) 



This in turn implies that the most general form of the hadronic tensor for unpolarized 
scattering is 



AttM 



M2 



(5.16) 



which defines the structure functions Wi and W2. 
Substituting Eq. ( p.l6| ) into Eq. ( |5.13| ) gives 



d^ 



a 



(2aE' 



dft'dE' \ Q2 



Wo cos 



2Wi sin^ 



(7 A/ 



W2 + 2Wi tan^ - 



(5.17) 



where 9 is the scattering angle of the electron in the lab system and ctm is the Mott cross 
section. 

We now calculate the hadronic tensor. If the spin of the target is in 1+3 dimensions 
this tensor is 



1 



{2S 



^P^^i2nmp, +P2-q- p)JlJa, 



(25 + 



5,si,S2 



(27r)64E;E^ 



47r(E 



(5.18) 



where S", Si, S2 are the spins, and the second line gives the result in the CM system. 
The volume integrals d^pi/ Ei and d^p2/E2 are covariant, insuring that the tensor is also 
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covariant. However, when we treat the hadronic degrees of freedom in 1+1 dimensions, 
consistency requires that the hadronic tensor also be evaluated in 1+1 dimensions. In this 
case Eq. ( p.l8| ) must be replaced by its 1+1 dimensional equivalent 



(2^+1) 



E 



{27r)HE[E2 
1 



ail 



(25+1) 



A\p[\iE[ + E^ 



(5.19) 



where the second line gives the result in the CM system with J7a^(±) corresponding to the 
cases when p'^ is parallel to q (+) or antiparallel (— ). Except for these two possibilities the 
delta functions completely fix the kinematics. 

To calculate the structure functions for the 1+1 dimensional models, we use Eqs. ( [4. 261) , 
and d+l). The definitions ( ^ToD and (|5l6|) are covariant, so the structure functions 
can be evaluated in any frame, and it is convenient to evaluate them in the CM frame of 
the outgoing hadrons. Furthermore, since q' and p' have components only in the t and z 
directions, gauge invariance insures that 



W 



1z <l'o 



(5.20) 



and that W^^ = W^^. Hence Wi and W2 can be extracted from 



p^OO _ y^zz 



AirM 



W2-W1 



(5.21) 



In the next two subsections we calculate the structure functions for the scalar and spin 1/2 
models in 1+1 dimensions. 

C. Structure functions for the scalar model 



Using Eqs. (|4.26| ) and ( |5.19| ), the structure functions for the scalar theory in the deep 
inelastic limit are 



ly^^ = 



2x Af^ 



where 



(2g-Pi-g')| 
(2g-pi-g')| 



{(I- Pi] 



{(I- Pi] 



{2q-pi -g2)2 



(Q ■ Pi)\ = {%E[ - q'zP'i)' = IQ' 



h — a 
b + a, 



I — X 



X 



(5.22) 



(5.23) 
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and 



a = 2p[q'^ 



2x 



X {ml + — 2M^x(l — x)) 



1 - 



2x 



2x 



x{l — 2x){m\ — mf) 

g2(i-x) 



(5.24) 



Hence 6 + a — Q'^/x, and the (— ) term vanishes in the deep inelastic hmit. The (+) term 
is finite however, giving 



Wi 

VW2 



X' 





X 

2^ 

, 

47r [miil — x) + m2X — M'^xil — x)Y 



X 



47r(l 
2q 



h — a 



X] 



(5.25) 



Wi vanishes in the deep inelastic limit, as expected for scalar particles. But vW2 scales 
in the deep inelastic region; its dependence on and u is replaced by dependence on the 
Bjorken scaling variable x alone. In the simple quark parton model, we expect 



Foix) 



xf{x) 



where f{x) is the probability of finding a quark with momentum fraction x 



dx f{x) 



1. 



Our simple model therefore predicts 



x(l — x) 



47r [mf (1 — x) + mix — M'^x{l — x)Y 



(5.26) 



(5.27) 



(5.28) 



which satisfies the normalization condition ( |5.271 ), as can be seen by comparing with 
Eq. 



D. Structure functions for the spin 1/2 model 

In the previous subsection we obtained the familiar result that the Callan-Gross relation 
is not satisfied by scalar quarks. This relation predicts that 



MW^iQ^u)^F,ix) = y{x), 



(5.29) 



where f{x) the same function which specified F2 in Eq. ( p.26| ). For scalar quarks the structure 
function Wi falls off with Q^. To obtain the relation (^.29|) , the quarks must to be treated 
as spin 1/2 fermions. [The diquarks, even when their orbital angular momentum is zero, 
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will have two spin states (0 and 1), and our choice of spin zero corresponds to the simplest 
possible wave function.] When carrying out the calculation in 1+1 D, we must assume that 
only the momenta are restricted to 1+1 D, but that the spin lives in the full 1+3 dimensional 
space. This assumption is fully consistent with the dynamical calculations of Sees. II and 
III. 

From the spin zero calculation we know that it is sufficient to evaluate the (+) matrix 
element only. In this case the incoming nucleon is traveling in the —z direction, and the 
final quark in the +z direction. We will calculate this matrix element in the CM system of 
the outgoing pair using helicity states, which for this case are 



u{pi, Ai) = ^Ep^+mi ^^^^^ j (5.30) 



where 



p = ~ 1 

^ Ep + M p 

Pi T mi 

Pi = — \ 1 

Ep^ + nil Pi 

With the help of these quantities, the x-component of the current ( p.21| ) becomes: 

M 



J: = ^{Ep, + mi){Ep + M) ^ 2Ajv [q\pi - p) + q\piP - 1)] 5a„,a, 



'^^ 2\n [Mx + mi] . (5.32) 



^/x[h- a] 



This gives 



^ AT^ {l-x){Mx + mif = ^ f( ^ r^QQ^ 

' S-kM [ml{l-x)+mlx- M^x{l-x)f 2M^^ ' ^ ' 



In a similar manner, it may be shown that W^*^ — W^^ in the deep inelastic limit, which 
is the Callan Gross relation [recall Eq. ( |5.21| )]. 

Note that the normalization condition for the spin 1/2 model given in Eq. ( |2.29|) confirms 



that this model is also correctly normalized. However, the /(x) obtained from the spin 1/2 
model does not go to zero as a; — > 0. 



E. Numerical results 



In Fig. |TU| the structure functions obtained from the scalar and spin 1/2 models are 
compared with the empirical momentum distribution xu{x) for u quarks in the proton. 
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FIG. 10. Quark momentum distributions: xu{x) with u{x) normalized to unity (dotted line), 
xf{x) for the scalar model with m'\/M'^ = 0.1 and mg/M^ = 2 (dashed line), and xfix) for the 
spin 1/2 model with m\/M'^ = 0.18 and rn^/M'^ = 2 (solid line). 

The empirical distribution was taken from Ref. and has been renormalized to unity for 
comparison. 

Note that the parameters in the momentum distributions predicted by the models can 
be chosen to give a qualitatively reasonable description of empirical results, but the models 
fail at both large and small x. The large x distributions fall like 1 — x, while the empirical 
distribution falls more rapidly (the fit gives (1 — x)^'^°^ compared to (1 — xY for the naive 
quark model). Similarly, the empirical distribution goes like at small x, while the 

spin 1/2 model goes like x and the scalar model goes like x^. Assuming the u and d 
quark distributions are identical, the resulting momentum carried by valence quarks is 3 x 
0.150 = 0.450 for the empirical distribution, 3 x 0.210 = 0.630 for the spin 1/2 model, and 
3 X 0.232 = 0.696 for the scalar model. In every case some momentum must be carried by 
glue or sea quarks, and the naive models presented here do not include these contributions. 
This will be investigated in a future work. 

VI. CONCLUSIONS 

In this paper we have shown how to extract the gauge invariant part of the pole term so 
that it gives the leading contribution to deep inelastic scattering. The gauge non-invariant 
part depends only on and is cancelled exactly by the rescattering term. While we have 
not shown it explicitly, the argument suggests that this procedure is quite general and works 
as follows. First add a (unique) term proportional to to the pole term which makes it 
gauge invariant. Then subtract the same term from all of the remaining terms (rescattering 
and interaction currents). It follows that the remaining terms are also gauge invariant, and 
we conjecture that the sum of all of these remaining terms will vanish in the deep inelastic 
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limit. Hence the modified pole term will give the exact result for DIS. 

More explicitly, imagine that the full current is the sum of a pole term Jp^jg and a 
remainder J/^m, and that it is gauge invariant. The first step is to write it as 

_ 7M I JfJ. 
pole ' '^rem 

= {jUe-^9') + i-Jrem + Cyqn , (6-1) 

where 

Note that a is uniquely determined (unless = 0, which cannot happen for DIS), and that 
the modified pole term 

■^modpole ~ "^polc ^2 ^ V^-'^) 

is gauge invariant by construction. We conjecture that the modified remainder terms (which 
are also gauge invariant) will vanish in the DIS limit, leaving ( |6.2[ ) to give the exact result 
for DIS. Note that our prescription ( |6.2|) is identical to the Landau prescription, and defined 
in Ref. §. 

The second result of this paper is the construction of toy models for the discription 
of DIS. These models give a qualitatively reasonable description of the phenomenology, as 
shown in Fig. |10|, but simplicity is their main virtue. We plan to use them to study many 
issues in the theory of DIS. 
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